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ABSTRACT 

Supcrmassivc black holes (BH) at the centres of galaxies can rapidly change their mass 
and spin by gas accretion and mergers. Using hydrodynamical cosmological simula- 
tions, with prescriptions for BH growth and feedback from Active Galactic Nuclei, we 
study the evolution of BH spins and we highlight the mechanisms responsible for driv- 
ing the amplitude and the direction of spins as a function of cosmic time. We find that 
in the high-redshift universe, galaxies maintain large values of gas accretion onto BHs. 
BHs tend to align their spins with the angular momentum of the surrounding gas and 
maximise their amplitude. On the opposite, at low redshift, as BHs get more massive, 
the contribution from binary coalescence to the total BH mass growth increases and 
tends to decrease the amplitude of spins and change their direction. 
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1 INTRODUCTION 

Supermassive black holes (BH) are commonly observed in 
the centres of galaxies. Their mass scales with the mass and 
the velocity dispersion of the host galaxy (Magorrian et al. 
1998; Tremaine et al. 2002; Haring & Rix 2004). This rela- 
tively tight correlation between the BH mass and its galaxy 
host properties has been interpreted as the signature of 
a self-regulating process operating from the BH onto the 
galaxy (Silk & Rees 1998; King 2003; Wyithe & Loeb 2003). 

It is now established that BHs can accrete gas from their 
surroundings and extract energy from the accretion process 
in the form of photons that couple to the optically thick gas 
to produce winds (quasar mode) or collimated jets (radio 
mode). Numerical implementations of the feedback of Ac- 
tive Galactic Nuclei (AGN) have demonstrated that the en- 
ergy released by BHs is able to shape the galaxy luminosity 
function in the bright-end (Croton et al. 2006; Bower et al. 
2006), and halt the cooling catastrophe operating in massive 
halos (Sijacki et al. 2007; Dubois et al. 2010; Teyssier et al. 
2011). 

The amount of radiative energy that is extracted from 
accretion events is closely related to the spin of the cen- 
tral BH the gas is accreting onto. Maximally spinning BHs 
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release more specific energy than non-rotating or counter- 
rotating BHs, which could have strong consequences for 
the self-regulated growth of BHs and the rate at which 
they accrete gas. Even though it is difficult to measure 
the spin of BHs from observations, several attempts have 
been made through X-ray spectra and produce contradic- 
tory results (Patrick et al. 2011; Brenneman et al. 2011; 
Reynolds et al. 2012). One recent measurement suggests 
that the BH with mass 2 x 10 6 M Q in the centre of NGC 
1365 is spinning at its maximum value (Risaliti et al. 2013). 

The direction of radio jets depends on the orientation 
of BH spins as jets emerge perpendicular to the inner part 
of the accretion disc which aligns with the spin. The capa- 
bility of jets to change directions is possibly a crucial phe- 
nomenon to help self-regulate the cooling flows in the centre 
of clusters of galaxies (Babul et al. 2012), even though ex- 
plosive jet outbursts could isotropise the energy output by 
shocks, sound waves or turbulence (Sanders & Fabian 2007; 
Dubois et al. 2010; Gaspari et al. 2011) 

Analytical or semi-analytical models based on pure 
dark matter simulations have tried to address the 
cosmological evolution of BH spins related to their 
host properties (Volonteri et al. 2005; Shapiro 2005; 
Volonteri et al. 2007; Lagos et al. 2009; Fanidakis et al. 
2011; Volonteri et al. 2012; Barausse 2012). Very few hydro- 
dynamical simulations of the spin evolution of BHs embed- 
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ded into their host galaxies have been performed: they are 
either performed in isolated contexts though with extremely 
high resolution (Maio et al. 2013) or neglect the role of gas 
accretion for spin evolution (Sijacki et al. 2009). Our aim is 
to explore further the spin evolution paradigm by evolving 
spins with gas accretion and coalescence that are dictated 
by their cosmological environment. 

In this paper, we explore the evolution of BH spins us- 
ing hydrodynamical cosmological simulations of galaxy evo- 
lution with BH growth and AGN feedback, i.e. where the 
properties of galaxies are treated self-consistently. BH spins 
are evolved through accretion of gas (tracking its angular 
momentum) and coalescences. In section 2, we describe the 
initial conditions and the physics (gas cooling, star forma- 
tion, stellar and AGN feedback) employed to follow the for- 
mation of galaxies in a ACDM universe. In section 3, we 
give the properties of our model for BH spin evolution by 
using properties extracted from hydrodynamical cosmologi- 
cal simulations. In section 4, we describe the predicted spin 
evolution from our model and how this quantity is driven 
by BH growth. We finally summarise and discuss our main 
results in section 5. 



2 NUMERICAL SET-UP 

2.1 Initial conditions 

We assume a ACDM cosmology with total matter density 
f^m = 0.272, baryon density fij = 0.045, dark energy den- 
sity Qa — 0.728, amplitude of the matter power spectrum 
as = 0.8 and Hubble constant Ho = 70.4 km s" 1 Mpc" 1 con- 
sistent with the WMAP 7-year data (Komatsu et al. 2011). 
The box size of our simulation is Lbox = 50/i _1 Mpc, with 
256 3 dark matter (DM) particles corresponding to a DM 
mass resolution of M rcs = 6.6 x 10 Mq. An extra zoomed- 
in simulation of a massive halo at z — 6 with 10 pc resolution 
is analysed further in Appendix A to test the robustness of 
sampling the gas angular momentum on kpc scales, along 
with a simulation of an isolated disc galaxy at similar reso- 
lution. 

The simulations are run with the Adaptive Mesh Refine- 
ment code RAMSES (Teyssier 2002). The evolution of the gas 
is followed using a second-order unsplit Godunov scheme for 
the Euler equations. The HLLC Riemann solver with a first- 
order MinMod Total Variation Diminishing scheme to re- 
construct the interpolated variables from their cell-centered 
values is used to compute fluxes at cell interfaces. Collision- 
less particles (DM, star and BH particles) are evolved using 
a particle-mesh solver with a Cloud-In-Cell interpolation. 
The coarse mesh (with minimum level of refinement equals 
8) is refined up to Aa; = 2.2 kpc (maximum level equals 15) 
using a quasi-Lagrangian strategy: when more than 8 DM 
particles lie in a cell, or if the baryon density is larger than 
8 times the initial DM resolution. The minimum cell size is 
kept roughly constant in physical size with redshift, i.e. an 
additional level of refinement is added every a oxp = n x 0.1 
(where n = 1, 2, 4, 8 and a exp is the expansion factor of the 



universe) up to level 15 at a a 



0.8. 



2.2 Physics of galaxy formation 

Gas is allowed to cool by H and He cooling with a con- 
tribution from metals using a Sutherland & Dopita (1993) 
model for temperatures above To = 10 4 K, which is the 
minimum temperature of the gas allowed through radiative 
losses. Heating from a uniform UV background takes place 
after redshift z rc i on = 10 following Haardt & Madau (1996). 
Metallicity is modelled as a passive variable for the gas ad- 
verted with the flow (whose composition is assumed to be 
solar) and is altered by the injection of gas ejecta during 
supernovae explosions and stellar mass losses. We assume a 
zero initial metallicity. The gas follows an adiabatic equa- 
tion of state for ideal monoatomic gas with adiabatic index 
7 = 5/3. 

The star formation process is modelled with a Schmidt 
law: p* = e*p/tff , where p* is the star formation rate den- 
sity, e* the constant star formation efficiency, and ts the 
local free-fall time of the gas. We choose a low star for- 
mation efficiency e* = 0.02 consistent with observations of 
giant molecular clouds (Krumholz & Tan 2007) and surface 
density relations of galaxies (Kennicutt 1998). Star forma- 
tion is allowed in regions exceeding a gas density thresh- 
old of po = 0.1Hcm~ 3 . The gas pressure is artificially en- 
hanced above p > po assuming a polytropic equation of 
state T = Tb(p/Po) K with polytropic index k = 4/3 to 
avoid excessive gas fragmentation. Feedback from stars is 
taken into account assuming a Salpeter initial mass function 
with a low-mass (high-mass) cut-off of 0.1 Mq (100 Mq), as 
described in detail in Kimm et al. (in prep.). Specifically, 
the mechanical energy from supernovae type II and stel- 
lar winds is taken from starburst99 (Leitherer et al. 1999, 
2010), and the frequency of supernovae type la explosions is 
computed following Greggio & Renzini (1983). The energy 
from SNe and stars is coupled to the gas with a kinetic im- 
plementation from Dubois & Teyssier (2008) until the age of 
stars reaches 50 Myr, where mass, momentum and kinetic 
energy profiles are imposed to mimic a Sedov blast wave. 
After 50 Myr, the energy is released thermally into the gas 
where the star is present. 



2.3 Model for BH growth and AGN feedback 

We use the same "canonical" AGN feedback modelling em- 
ployed in Dubois et al. (2012). BHs are created at loci where 
the gas density is larger than the density threshold for star 
formation po with an initial seed mass of 10 5 Mq. In order 
to avoid the formation of multiple BHs in the same galaxy, 
BHs are not allowed to form at distances smaller than 50 
kpc from any other BH particle. The accretion rate onto 
BHs follows the Bondi-Hoyle-Lyttleton (Bondi 1952) rate 
A/bh = 47raG 2 Afl H p/(c? + u 2 ) 3/2 , where M BH is the BH 
mass, p is the average gas density, c s is the average sound 
speed, u is the average gas velocity relative to the BH veloc- 
ity, and a is a dimensionless boost factor with a = (p/po) 2 
when p > po and a = 1 otherwise (Booth & Schaye 2009) in 
order to account for our inability to capture the colder and 
higher density regions of the ISM. The effective accretion 
rate onto BHs is capped at the Eddington accretion rate: 
Mvdd — 47rGMBHWip/(e r aTc), where ot is the Thompson 
cross-section, c is the speed of light, m p is the proton mass, 
and e r is the radiative efficiency, assumed to be equal to 
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e r = 0.1 for the Shakura & Sunyaev (1973) accretion onto a 
Schwarzschild BH. Note that the radiative efficiency of ac- 
cretion is a function of the spin parameter of the BH. Here, 
we assume that value to be fixed as the BH spins are not 
followed on-the-fiy, but are a post-processed quantity from 
cosmological runs. 

In order to avoid spurious oscillations of the BH in the 
gravitational potential well due to external perturbations 
and finite resolution effects, we introduce a drag force that 
mimics the dynamical friction exerted by the gas onto a 
massive particle. This dynamical friction is proportional to 
Fdf = / g as47rap(GMBH/c s ) 2 , where / gas is a fudge factor 
whose value is between and 2 and is a function of the 
mach number M = u/c s < 1 (Ostriker 1999; Chapon et al. 
2013), and where we introduce the boost factor a for the 
same reasons than stated above. 

The AGN feedback is a combination of two differ- 
ent modes, the so-called radio mode operating when \ = 
Mmi/M-Ed& < 0.01 and the quasar mode active otherwise. 
The quasar mode corresponds to an isotropic injection of 
thermal energy into the gas within a sphere of radius Aa;, 
at an energy deposition rate: -Bagn = £{£tMbhc 2 , where 
£f = 0.15 for the quasar mode is a free parameter cho- 
sen to reproduce the Mbb-M^, Mbh-oIj, and BH density in 
our local Universe (see Dubois et al. 2012). At low accretion 
rates on the other hand, the radio mode deposits the AGN 
feedback energy into a bipolar outflow with a jet velocity 
of 10 4 km.s -1 into a cylinder with a cross-section of radius 
Aa; and height 2 Aa; following Omnia et al. (2004) (more de- 
tails about the jet implementation are given in Dubois et al. 
2010). The efficiency of the radio mode is larger with a — 1. 



3 MODEL OF BH SPIN EVOLUTION 

3.1 Gas accretion 

BHs are allowed to be spun up by accretion of gas through 
the following expression (Bardeen 1970): 
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where M ra tio = M^ /Mg H , the n superscript stands for the 
values at time t n , and i?i sco is the radius of the innermost 
stable circular orbit (ISCO) defined as (in reduced units): 



Risco/R s = 3 + Z 2 ±[(3- Zi)(3+ Z1 + 2Z 2 ) 
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(2) 



where the gravitational radius i? g is defined as half of the 
Schwarzschild radius of the BH, -Rbh, and Z\ and Z2 are: 
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For positive spins, co-rotating with their accretion disc, 1 ^ 
r'isco < 6, while for negative spins, counter-rotating with 
their accretion disc, 6 < risco ^ 9, and ri sco = 6 for non- 
spinning BHs. 

The equation that governs the evolution of BH spin 
through direct accretion of gas is true when BH spin and disc 
angular momentum are perfectly aligned (or anti-aligned), 
but in the most-general case, misalignment occurs. In the 
misaligned case, the accretion disc experiences a torque due 



to the Lense-Thirring effect that causes the accretion disc 
to precess about the spin axis of the BH. For large enough 
viscosities the innermost parts of the disc are forced to ro- 
tate within the equatorial plane of the disc and a warped 
disc is created. We can define the total angular momen- 
tum of the system {disc+BH} as Jtot = Jd + Jbh. The 
values of the angle 8 between Jbh and Jd are between 
— 1 ^ cos# ^ 1 where the two-extrema correspond to anti- 
aligned and aligned cases respectively. 

The result of the Lense-Thirring precession is that the 
BH and disc angular momentum ends up being aligned 
or anti-aligned with the total angular momentum. The 
case for anti-alignment of the BH with the disc requires 
that (King et al. 2005) 

COS 0<_ JfL_, (5 ) 

^Jbh 

thus, for cos# ^ 0, they always align, while for cos 9 < 0, 
they eventually anti-align if the ratio Jd/ Jbh is sufficiently 
small compared to cos#. 

The difficulty is to access to the magnitude of Jd as we 
cannot resolve the accretion disc. We therefore assume the 
tilted solution for a thin accretion disc (Shakura & Sunyaev 
1973; Scheuer & Feiler 1996; Natarajan & Pringle 1998; 
Perego et al. 2009), characterized by a viscosity v\ = 
QiCsfGMBH/r 3 ) 1 , where a t is a parameter < 1, and c s 
is the sound speed. In a misaligned disc under the effect 
of Lense-Thirring precession a natural scale is the warp ra- 
dius, that marks the transition between an equatorial in- 
ner disc and a misaligned outer disc: the direction of the 
angular momentum of the inflowing material changes direc- 
tion as it passes through the warp. Only material within the 
warp radius can effectively transfer its angular momentum 
to the BH (Volonteri et al. 2007). This is the relevant scale 
to estimate whether alignment or ant i- alignment occur. In 
a Shakura & Sunyaev (1973) thin accretion disc, one can 
write the ratio of the warp to Schwarzschild radius as: 

-5/8 
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where £ r ,oi is the radiative efficiency normalised to 0.1, l>\ 
and V2 are the kinematic viscosities horizontal and perpen- 
dicular to the equatorial plane of the disc: vi is the viscos- 
ity responsible the warp propagation, while v\ is the vis- 
cosity responsible for driving accretion and transferring an- 
gular momentum. We choose a value for at, 01 = at/0.1 = 
1 (King et al. 2007), and the ratio of v 2 /vi = 2(l+7a t )/(4+ 
a? )/ Q t is of the form given by Ogilvie (1999) (equation 145). 
With these choices, (^2/^1) ~ 85 (note that for small at, 
v 2 ~ ui/ott) and: 
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The BH angular momentum is simply Jbh = 
3/2 1/2 
aMgJj R B b i a 11 *! ^ ne disc angular momentum can now be ex- 
pressed as Jd ~ Md(-Rwar P )A/ B (j i?warp. The disc mass within 
-Rwarp is Md(-Rwarp) = Mt l/1 (R vra . Tp ) where t vi is the viscous 



timescale for radial propagation: 
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The ratio of disc to BH angular momentum can, then, 
be written as: 
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Equation (9) allows one to evaluate if anti-alignment occurs, 
and provides a value for the magnitude of Jd- We refer the 
reader to Dotti et al. (2013) for a detailed discussion of how 
the ratio Ji/Jbh determines the behaviour of alignment and 
spin evolution. 

For the direction of Jd, we make the assumption that 
the accretion disc is aligned with the gas angular momen- 
tum, J g , measured in the surroundings of the BHs as ex- 
tracted directly in the simulation (but their magnitude is 
not equal). Note that it is potentially a strong assumption 
because, if the gas is sufficiently turbulent, the inner dense 
cold (unresolved) molecular gas surrounding the BH can 
misalign with the large-scale galactic disc. In Appendices A 
and B we provide additional details and validate our tech- 
nique and assumption using high-resolution simulations that 
purposely study the degree of consistency between the an- 
gular momentum on large scales and small scales, as well as 
the dependability of this assumption. 

Thus, we update the orientation of the BH angular mo- 
mentum after one accretion event (in order to avoid over- 
estimating Jd we divide each simulation coarse time step, 
~ 2 Myr, in a sequence of At = t vi ) with J^ 1 = Jtot = 
•^BH + ^dJg f° r the alignment case, where the small j stands 
for angular momentum unit vectors, and where Jd is pro- 
vided by equation (9). If the criterion for anti-alignment is 
met, J^ 1 becomes —Jtot if Jtot- Jbh > 0, this occurs when 
the total angular momentum dominated by the BH (case b in 
Fig. 1 in King et al. 2005). Otherwise, if the anti-alignment 
condition is satisfied and Jtot -Jbh < 0, then the total angu- 
lar momentum is dominated by the disc and it is anti-aligned 
with respect to that of the BH, therefore JJ^h 1 = Jtot (case 
d in Fig. 1 in King et al. 2005). 



3.2 Coalescence of BHs 

We track the change of the spin during the coalescence of two 
BHs using the analytical fit of Rezzolla et al. (2008) derived 
from relativistic numerical simulations of BH binaries. The 
final value of spin of the BH remnant is 
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where ai is the spin vector of the most massive BH pro- 
genitor, <i2 is the spin vector of the least massive BH pro- 
genitor, q — M2/M1 ^ 1 the mass ratio of the binary, and 
£ = lj(M\M-2), with I the difference between the orbital an- 
gular momentum L when the two BHs are widely separated 
and the angular momentum J gw extracted from gravita- 
tional waves before the final coalescence I — L — J ew - 



Rezzolla et al. (2008) provides a simple analytical ex- 
pression for the norm of I 
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where <j>\ (<f> 2 ) is the angle between a\ (02) and £, 
and n — 5/(1 + q) 2 ■ We made the same assumption 
than Rezzolla et al. (2008) to get the direction of £ by 
imposing collinearity with the orbital angular momentum 
L — L\ + Z/2 before coalescence, which is a vector directly 
measured in our simulations. To get access to L, we com- 
pute the centre of mass of the binary r com and we evaluate 
the angular momentum of the two BHs relative to the cen- 
tre of mass with L\ = Mi(n — r com ) x (vi — » mm ). The 
numerical constants involved in equation (11) are equal to 
s 4 = -0.129, s 5 = -0.384, to = -2.686, t 2 = -3.454, and 
£3 = 2.353. 



3.3 Updating spins 

The whole evolution of BH spins is post-processed over the 
numerical simulation results. For this purpose, we print for 
each coarse time step (synchronisation of all levels at the 
minimum level of refinement) of the simulation the proper- 
ties of BH particles that allows to derive the evolution of 
their spin. We assume that the initial spin of BHs is zero 
at creation. We further assume that the thin disc solution 
for spin evolution applies at any accretion rate over Edding- 
ton x- It seems to be a fairly good approximation as BH 
spins mostly change their value by accretion during strong 
accretion events (of the order of \ < 1). The maximum spin 
value that we authorise is a max = 0.998, due to the capture 
by the central BH of radiated photons emitted from the ac- 
cretion disc (Thorne 1974). Note that a max = 0.9 is also 
possible for thick discs taking into accounts the magneto- 
hydrodynamical effects on the transfer of gas angular mo- 
mentum (Gammie et al. 2004), while Sadowski et al. (2011) 
suggest a ma x = 0.9994 if BH accretes gas at super-Eddington 
within a low-viscosity a t = 0.01 accretion disc. 

We insist on the fact that the BH spin is a post- 
processed quantity from the hydrodynamical cosmological 
simulations. For this reason, we neglect the effect of grav- 
itational recoil on BH particles during the coalescence of 
BH binaries, which kicks the remnant BH depending on 
the amplitude and configuration of their spins before coa- 
lescence (Koppitz et al. 2007). In a future work, we will im- 
plement the scheme described above directly in the RAMSES 
code, which will make possible to determine the efficiency 
of AGN feedback and the direction of jets from the spin 
amplitude and direction. 

Note that a range of theoretical models assume that 
gas accretes onto BHs in random directions due to the ca- 
pacity of the cold gas of the interstellar medium to drive 
turbulence (Hopkins et al. 2012) or due to the formation of 
a self-gravitating accretion disc around BHs (e.g. King et al. 
2008). Therefore, we run extra models where some level of 
randomness in the direction of the accreted gas angular mo- 
mentum is added explicitly to the scheme, either continu- 
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Figure 1. Example of the evolution of different properties of a BH of mass 5 X 10 9 Mq and which mass has been built-up by 73 per 
cent through mergers at z = 0. From top to bottom: BH mass evolution, BH spin evolution (red points indicate coalescence with its 
BH mass ratio), cosine of the angle between BH spin and gas angular momentum 8, z-cartesian component of the normalised BH spin 
vector, z-cartcsian component of the normalised gas angular momentum (both calculated with respect to the fixed reference frame of the 
simulation box), accretion rate over Eddington \- 



ously or activated only during short starburst triggered dur- 
ing galaxy mergers. The results of the models with random 
gas directions are discussed in Appendix B. 



4 RESULTS 

4.1 General properties of the gas flow near BHs 
from the simulations 

Before discussing in detail the results of the spin models, we 
briefly describe the properties of the gas flows near BHs and 
how they come to influence the angular momentum of the 
BHs. These results are derived directly from the simulations, 



before any post-processing is applied. In Fig. 1, 2 and 3 we 
show some individual BH histories that highlight the main 
findings. In this section we focus only on the two bottom 
panels, that show the degree of temporal coherence of the 
gas angular momentum through j gas , z , which is the cartesian 
component along the (fixed) z-axis of the simulation box, of 
the gas angular momentum unit vector, and the accretion 
rate in Eddington units, \- 

Fig. 1 shows the evolution of a BH hosted in what 
becomes by z = a central galaxy of a cluster of galax- 
ies. The BH mass is 5 x 10 Mq and lies in a host galaxy 
with a stellar mass of 9.7 x 10 11 Mq, and host halo mass of 
Mvir = 1-6 x 10 14 M . At high redshift, z > 2, the BH regu- 
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Figure 2. Same as Fig. 1 zoomed on the redshift range z = [0,0.2], 



larly accretes gas at or close to its Eddington limit. The gas 
possesses a high degree of coherency, with small variations 
of the degree of alignment with the galaxy total angular mo- 
mentum. At low redshift, z < 2, the galaxy has consumed 
most of its gas through star formation, and blown some gas 
away through feedback (mostly due to AGN feedback in 
such a massive elliptical). Fig 2 highlights this new phase, 
by zooming in on < z < 0.2 for the same BH. The only 
available gas is hot plasma, at temperatures of a few 10 7 K 
and densities of 0.01-0.1 cm" 3 , that does not settle into a 
well structured disc. The rapid variations of the gas angu- 
lar momentum direction are driven by the local turbulence 
and bear no connection to the rotation axis (if any) of the 
galaxy. Accretion rates at low redshift can vary by several 
decades from a few x = 10~ 6 over long period of time up to 
X = 1 with short burst durations. In general, however, the 



mean accretion rates are very low (except for brief bursts due 
small, short-lived cooling flows). The tiny accretion rates are 
dictated by the absence of any cold star-forming gas in the 
galaxy. Gas accretion, in this phase, does not contribute sig- 
nificantly to the BH evolution, nor it mass, neither its spin, 
as detailed below. 

In Fig. 3, we show instead the evolution of a massive 
disc galaxy of M s = 5 x 10 10 M Q in a Af vir = 2 x 10 12 M Q 
halo hosting a Mbh = 10 M© BH at z = 0. This galaxy lies 
in the field (its host halo is not a sub-halo of a more mas- 
sive halo) and contains in permanence a reservoir of cold 
star-forming gas supplying accretion for its central BH (and 
for star formation). In this case we see that the gas main- 
tains a coherent angular momentum down to low-redshift, 
while at the same time the reservoir of cold gas keeps on 
feeding the BH at relatively high-levels. The bottom panel 
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Figure 3. Same as Fig. 1 for a BH of mass 10 8 Mq and which mass has been built-up by 6 per cent through mergers at z = 0. 



shows a general negative gradient of the accretion rate with 
decreasing redshift (from x~0.1atz = 3 — 4 to x ~ 10 -3 
at z < 1), highlighting that gas is still being consumed over 
cosmic time. Overall, it is gas accretion that drives this BH's 
evolution, both its mass, and its spin. 

These properties of the galaxy gas can help us under- 
stand how its angular momentum, temperature and density 
affect (or in some case do not affect) BH spins. To assess 
the robustness of these results we describe in Appendix A 
two simulations at much higher resolution (fO pc), one of a 
high-redshift halo in a cosmological simulation, the other of 
an isolated halo. The former allows us to estimate the effects 
of galaxy mergers and of cosmic feeding, the latter allows us 
to highlight the specific influence of AGN feedback and of 
the dynamics of gas clumps, all this down within nuclear 
scales. In all cases we find that the gas angular momentum 



at scales of tens of parsec shows a significant degree of cor- 
relation with the gas angular momentum on galactic scales, 
thus validating our assumptions. 



4.2 Examples of BH mass and spin evolution 

In order to understand the connection between the BH mass 
growth and spin, we inspect the redshift evolution of the 
two BHs described in section 4.1. Fig. 1 shows one of the 
most massive BHs in our simulation, hosted at z = in the 
central galaxy of a cluster of galaxies. Fig. 1 shows the BH 
mass as a function of redshift with its related spin evolution, 
angle between BH spin and gas angular momentum, the z- 
cartesian component of the BH and gas normalised angular 
momenta, and the BH accretion rate in Eddington units. At 
high redshift, z > 2, the BH steadily accretes gas close to 
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its Eddington limit, with strong coherence of the accreted 
gas angular momentum, because the direction of the gas an- 
gular momentum with redshift does not vary much. In this 
situation, the large values of \ force the BH to re-align with 
its surroundings (cos 8 ~ 1) and the BH spins at its maxi- 
mum a ~ 1. Mergers between BHs (and their host galaxies) 
are present in the high redshift evolution, and, if BH mass 
ratios are sufficiently large, they change the value of the BH 
spin (e.g. at z = 2.9 from a ~ 1 before merger to a ~ 0.8 
after merger). However, as mergers are still extremely gas 
rich, a rapid accretion of cold gas perdures after the coales- 
cence of the two BHs, and the spin value gets soon back to 
its maximum a — 0.998. 

The mechanism at low redshift (z < 2) is orthogonal to 
the picture describing the spin evolution at high redshift. 
As gas is depleted through star formation and feedback, 
spin decouples from the surrounding gas angular momen- 
tum. The angle between the two angular momentum direc- 
tions (gas and BH) strongly varies with redshift, essentially 
due to rapid variations of the gas angular momentum direc- 
tion driven by the local turbulence of the hot plasma (see 
Fig. 2). Conversely, the BH spin direction is much more sta- 
ble with long periods of time without changing direction 
(several hundreds of Myr) because the rates of gas accretion 
onto the BHs are too small to realign the BH spin with the 
local gas angular momentum. Since for high-mass BHs the 
total angular momentum of the BH+disc system, Jtot, is 
dominated by the BH's spin, therefore, although the direc- 
tion of Jd changes erratically, the re-orientation of the BH 
is limited, as its spin direction is already close to the direc- 
tion of Jtot (cf. Dotti et al. 2013). The rapid changes in spin 
magnitude and directions are due either to BH coalescences 
that redistribute the final BH spin direction (e.g. at z = 0.2 
where a merger with BH mass ratio of 1:2 happens), or to 
the formation of a small cooling flow of gas onto the central 
galaxy (e.g. at z = 1 where % ^ 10 -5 but no mergers are 
detected) . 

In Fig. 3 we show the results for a A/bh = 10 8 M BH 
within a massive disc galaxy. This field galaxy experiences 
very few mergers (the fraction of BH mass gained by coa- 
lescences is 6 per cent) and remains rich in cold gas down 
to 2 = 0. The difference with the previous BH is striking, 
the BH maintains a spin parameter close to its maximum 
a — 0.998 because of the remarkable absence of significant 
mergers (except at z — 3.4) and the moderate levels of ac- 
cretion (x — 10~ 3 -10 -2 ) maintained over a Hubble time in 
a coherent fashion (direction of the gas angular momentum 
varies little with time). 

These two BHs are extreme cases with one essentially 
built-up by mergers and the other one built-up by accretion. 
Indeed, it exists a range of BHs with moderate fractions 
of mass gained through mergers for different values of BH 
masses. In the following, we will pinpoint how the mass built- 
up of BHs determines the spin of these compact objects. 



4.3 BH spins and their relation to gas accretion 
and mergers 

Fig. 4 shows the spin value as a function of the BH mass 
at z = 0. Note that BHs with masses below 10 M© are not 
represented, because we consider that they are too close to 
the chosen initial seed BH mass, and they are still affected 
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Figure 4. Top: Average fraction of BHs with negative spins as a 
function of the BH mass; bottom: Median of the absolute value of 
the spin as a function of the BH mass (solid line), with 20 and 80 
percentile values (dashed lines). Measurements are made at z = 0. 
The red curve corresponds to the median of the absolute value of 
the spin for a model of spin evolution where BH coalescence does 
not modify the spin of the remnant. 



I.U 


■ 


0.8 


- 




. 


- 


0.6 






- 










. 








/ 




■ 


0.4 












/ [ : 


0.2 
















- 












. 


00 


-^--M-"" 









10 7 



10" 

Mbh (M sun ) 



10" 



Figure 5. Average fraction of BH mass gained through coales- 
cence with companions as a function of the BH mass at z = 0. Er- 
ror bars correspond to the standard dispersion. Contribution from 
BH coalescence to the final mass is more important for the most 
massive BHs. Lower mass BHs gain mass preferentially through 
gas accretion from their environment. 



by the mass resolution sampling. The least massive BHs we 
consider (10 6 M© < Mbh < 1O 7 M0) have large spins, close 
to \a\ - 0.998, while the most massive BHs (M BH ^ 10 s M Q ) 
tend to lower spin values, with \a\ ~ 0.7 — 0.8 for BHs with 
mass Mbh — 10 9 Mq . Coherent gas accretion onto BHs, i.e. 
gas angular momentum aligned with BH spin, is the rea- 
son for the spin values approaching a = 0.998: gas in disc 
galaxies is set-up on almost steady orbits, and contributes 
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Figure 6. Top: Average value of the Eddington ratio as a func- 
tion of redshift for selected BHs with masses at z = above 
M BH > 1O 8 M (black), A/ BH > 10 9 M Q (red), and 10 8 sj 
A^BH ^ 10 9 Mq (blue) with their standard dispersion (error 
bars). Bottom: Average value of the fraction of BH mass gained 
through mergers. The average accretion rate onto BHs is steadily 
decreasing over time, and at redshift z ~ 2 the most massive 
BHs switches from a quasar mode (x ^ 0.01) to a radio mode 
(x < 0.01). The contribution of mergers to the BH mass is in- 
creasing with time (decreasing redshift), starts to be significant 
by z = 2-3, and the most massive BHs have the largest mass 
fraction of mergers. 



constructively to the spin amplification. Note that a non- 
negligible fraction of BHs, between 10 and 20 per cent, have 
a negative spin value that also tend to a — —0.998 for low 
mass BHs and to a more intermediate value a ~ —0.7 for 
the most massive BHs. Midsize mass BHs (Mbh < 10 8 Mq) 
with negative spins have the gas angular momentum vec- 
tor flipped with respect to the spin direction because feed- 
back from the central AGN removes gas from the central 
region, while at the same time their typical levels of accre- 
tion are sufficiently low that the magnitude of the spins does 
not change to intermediate spin negative values through ac- 
cretion of counter-rotating material. At larger masses, spin 
amplitudes deviate significantly from unity and this may be 
due to either: (i) gas angular momentum decoherence with 
respect to the BH spin or (ii) coalescence of BHs. 

Coalescences of BHs can rapidly change the BH spin 
amplitude and direction. We measure the fraction of BH 
mass that has been gained through coalescences with com- 
panions (as opposed to growth by direct gas accretion) 
/merge = AM m erge (< t)/MsH, where AM me rge (< t) is the 




Figure 7. Same as Fig. 7 for the average spin parameter. Spins 
are maximum at high redshift when accretion drives the growth 
of BHs and the spins decrease at low redshift when binary coa- 
lescence start to contribute significantly to the final BH mass. 
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Figure 8. Same as Fig. 7 for the average radiative efficiency. Ra- 
diative efficiencies follow the evolution of the spins with redshift, 
with maximum efficiencies at high redshift and reduced efficien- 
cies at low redshift because of lower spin values. 



total integrated mass gained by mergers for a given BH, 
as a function of BH mass at z — in Fig. 5. Note that 
this definition does not take into account the fact that the 
merged companion of a given BH has grown by direct gas 
accretion, all its mass before coalescence is added to the 
merged mass M mer ge- Low mass BHs Mbh < 10 8 Mq have 
gained most of their mass through direct accretion of gas, 
/merge < 0.2, while the most massive BHs are dominated 
by mergers, /merge > 0.5 for BHs with Mbh > 10 9 Mq (see 
also Malbon et al. 2007; Volonteri & Ciotti 2012). It appears 
that the most massive BHs at z = that have grown mostly 
by coalescences have lower spin values as already seen in 
Fig. 4, and thus, galaxy mergers at low redshift that carry 
little gas, are responsible for spinning down BHs. We tested 
this scenario by deactivating the spin evolution during BH 
coalescence, i.e. BHs mass is still updated but not their spin, 
and we see in Fig. 4, that the red curve is extremely close to 
a — 0.998 even for the most massive BHs. It clearly demon- 
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strates that BH mergers are responsible for spinning down 
BHs at low redshift. 

Fig. 6 shows the evolution of the Eddington ratio with 
redshift for a sub-sample of BHs 1 . We selected BHs with 
masses above 10 8 Mq at z — and followed their progen- 
itors at higher redshift. The progenitors of these massive 
BHs accrete gas close to the maximum Eddington rate at 
high redshifts, and release their energy into quasar mode. 
At z ~ 2 these BHs switch from a quasar mode to a radio 
mode of AGN feedback (\ < 0.01), and the Eddington ratio 



reaches x 



10" 



10 at z = 0, which correspond to a 



typical growth rate of tEdd/x ~ 1000 Gyr, where iEdd = 45 
Myr is the e-folding time of a BH accreting at the Edding- 
ton rate (assuming e r = 0.1). It means that, at low redshift, 
the most massive BHs stop growing by accretion of gas from 
their surroundings, they can eventually get more massive by 
mergers. Indeed, the fraction of mass gained through merg- 
ers is continuously increasing with time, and represents 30 
per cent of the mass built-up of BHs at z = for those with 
Mbh > 10 8 Mg, as we discussed above (cf. Fig 5). 

An interesting point is that more massive is the BH the 
faster the quenching of gas accretion (relative to Eddington) 
is reached, e.g. progenitors of 10 9 Mq BHs at z = have 
lower Eddington ratio than progenitors of 10 8 Mq BHs. The 
same is true for the fraction of mass gained through mergers: 
mergers start to be significant (/merge > 0.1) earlier (at z ~ 3 
for Mbh > 10 9 M instead of z ~ 2 for 10 8 M < M BH < 
1O 9 M0), and represent a more significant fraction of the 
mass growth at z — (/merge — 0.55 instead of 0.3). 

We compute the evolution of the average spin parameter 
as a function of redshift for the same sample of supermassive 
BHs (see Fig. 7). At high redshift z > 2, the progenitors of 
the most massive BHs at z — have spin close to their max- 
imum value a — 0.998. At lower redshift (below z < 1 — 2), 
the value of the spin decreases with time as BHs experience 
more mergers that change spin directions and spin down 
BHs. As gas accretion onto BHs is quenched efficiently at 
low redshift, gas cannot spin back up the BHs efficiently. 

From the values of spin, it is possible to evaluate the 
radiative efficiency (which is kept as a constant in our feed- 
back model and equal to e r = 0.1) of the accretion process 
which writes e r = 1 — Ei BCO = \fl — 2/(3n SC o), where Ei sco 
is the specific energy of a gas particle in the innermost sta- 
ble circular orbit. Thus, for a maximum spin of a = 0.998, 
the radiative efficiency is equal to e r = 0.321 (e r = 0.057 
for a = and e r = 0.038 for a = -0.998), which is the 
exact value reached for BHs at high redshift when they are 
maximally spinning (see Fig. 8). While spins are decreasing 
with redshift, the radiative efficiencies are also decreasing, 
reaching e r = 0.11 for Mbh > 10 9 Mg at z — and larger 
values for less massive BHs. 

The high values of spin, and, thus, of radiative efficien- 
cies measured at high redshift have possible important con- 
sequences for the growth of high-redshfit BHs, as the most 
massive ones need to grow by several orders of magnitude 
to reach 10 9 Mq at z — 6 in less than a billion year to 



A glitch in the Eddington ratio is visible at z = 4 and is due to 
one extra level of refinement triggered at that moment, increasing 
the force sampling in galaxies and, thus, the density of the cold 
gas, which in turns produce larger values of the gas accretion rate. 
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Figure 9. Number-weighted distributions of angles between BH 
spins and galactic angular momentum (defined by their stellar 
content) as a function of the BH mass at z = 0. The solid line 
is the average of the distribution. Midsize mass BHs have spin 
aligned with their host galaxy, while the most massive BHs arc 
randomly oriented. 



explain observed powerful quasars in the high redshift Uni- 
verse. Given the timescale, they need to accrete gas at the 
Eddington rate, but given an efficiency of e r = 0.321 the 
growth time scale is tEdd(e r = 0.321) = 145 Myr (to be com- 
pared to 45 Myr for e r = 0.1), and, thus would require a 
seed mass BH of 10 6 Mq possibly somewhat unrealistic. If 
the spin is so large even at high redshift, this would favour 
a super-Eddington mode of accretion, or that the efficiency 
in producing radiation is decreased by, e.g., jet production 
(Ghisellini et al. 2013). 

Some empirical constraints from Li et al. (2012) 
and Shankar et al. (2013), using continuity equation ap- 
proaches, suggest that the average radiative efficiency de- 
creases with decreasing redshift (in agreement with these re- 
sults). However, Shankar et al. (2013) favours lower typical 
radiative efficiencies, in general, between 0.05 and 0.14 (see 
also Merloni & Heinz 2008). This models also prefer an in- 
creasing radiative efficiency with increasing BH mass (op- 
posite to our results). On the other hand, Li et al. (2012) 
find almost no dependence of the radiative efficiency with 
BH mass at z < 1 and an increasing radiative efficiency 
with mass at z > 1 (see their Fig. 6). Regarding the red- 
shift dependence, their results for the most massive BHs 
(logA^BH > 8.5) is in perfect agreement with ours, while at 
lower masses (7.5 < log Mbh < 8.5) the result is opposite to 
ours: no redshift dependence, and low values of the radiative 
efficiency, e r ~ 0.05. 



4.4 Spin amplitude and orientation with their 
host galaxy 

We saw that mergers are responsible for spinning down BHs, 
and that the most massive BHs with larger contribution 
from mergers to their mass are those with lower spins. Fig. 9 
shows the distribution of angles between BH spins and the 
stellar angular momentum ^BH.gai of their host galaxy as a 
function of BH mass. Midsize nrass BHs with 10 7 < Mbh < 
10 8 Mq show a high level of alignment with their host galaxy, 
as gas accretion is mostly responsible for both BH and stellar 
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Figure 10. Number-weighted distributions of angles between BH 
spins and galactic angular momentum (defined by their stellar 
content) as a function of the gas fraction in their host galaxy at 
2 = 2 (top), z = 1 (middle), z = (bottom). Solid lines are the 
average of the distribution. Gas poor galaxies host BHs whose 
spin is misaligned with the angular momentum of the stars (but 
not totally random), while gas-rich galaxies host strongly aligned 
BHs. 



mass growth in this mass range. As BHs get more massive, 
spins are more randomly oriented with respect to their host 
galaxy angular momentum and tend towards cos #BH,gai = 0. 
Mergers can rapidly change the orientation of galaxies, as 
do coalescences for BHs but they do not necessarily follow 
the same evolution while they merge. BHs radiate a part 
of their orbital angular momentum in gravitational waves 
that change the final orientation of the BH remnant relative 
the galaxy angular momentum. Galaxies can also exchange 
angular momentum with their environment (gas and DM) 
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Figure 11. Number-weighted distribution of BH mass fractions 
gained through coalescence (top) and distribution of BH spins 
(bottom) as a function of the gas fraction in their host galaxy 
at z = 0. Solid lines are the average of the distribution. BHs 
whose mass is mostly gained through coalescence instead of direct 
gas accretion are preferentially hosted in gas-poor galaxies, which 
results in BHs with lower spin values. 



by long-distance gravitational interactions. Low mass BHs 
show a significant deviation from cos^BH.gai = 1, which is to 
be attributed to resolution effects. A BH with mass Mbh = 
10 6 M Q lie typically in host galaxy of M s = 10 9 M which is 
composed of ~ 50 stellar particles (while the hosts of BHs 
with mass Mbh = 10 7 Mq are composed of ~ 500 particles, 
and so forth). 

We investigate if there is any correlation between the 
amount of gas in galaxies and the spin amplitude or ori- 
entation of their BHs. We define the ratio of gas in the 
galaxy by Mism/(Mism + M a ), where M a is the total stellar 
mass of the galaxy obtained with a halo finder (AdaptaHOP, 
see Aubert et al. 2004; Tweed et al. 2009) performed on star 
particles, and Mism is the star forming gas within the galac- 
tic radius with gas densities n > 0.1 cm~ . Fig. 10 shows the 
distribution of angles between BH spins and their stellar an- 
gular momentum #bh, g ai of their host galaxy as a function 
of the gas fraction in the galaxy for BHs with masses above 
Mbh > 10 7 M Q . There is some correlation of 0bh,gal with 
the gas content in the galaxy: BHs in gas-rich galaxies have 
spins aligned with the galaxy angular momentum, and BHs 
in gas-poor galaxies have more randomly orientated spins 
though with some degree of correlation, cos#bh.gai = 0.5 for 
the galaxies with zero star-forming gas. Gas-rich galaxies 
provide both strong accretion onto BHs and actively star- 



12 Y. Dubois et al. 



forming regions through the same gas, thus contributing 
to aligning together BHs, stars and gas. Viceversa, galax- 
ies with low gas content are the systems less effective at 
realigning spins with stellar angular momentum, and these 
galaxies are those which host BHs whose fraction of mass 
gained through coalescence is the largest (see Fig. 11, top 
panel). However the orientation of spins is not purely ran- 
dom even for the gas-poor galaxies. In the case of the most 
massive galaxies, they show some degree of alignment be- 
cause they experience intense mergers that efficiently reori- 
entate spin and galactic angular momentum. Lower mass 
galaxies, instead, have fewer mergers and keep the BH spin 
aligned with their stellar angular momentum as they both, 
the BH and the stellar component, grow from the same gas. 
Galaxies with low gas fractions exhibit lower degree of align- 
ment with their central supermassive BHs at higher redshift. 
This is essentially an effect of the quiescent depletion of gas: 
gas-rich galaxies at high redshift, that host strongly aligned 
BHs, consume their gas through the star formation process 
(and partially because of AGN feedback, see Dubois et al. 
2012) and end up with lower gas fractions at low redshift, 
with BHs still aligned to their hosts if they have not endured 
a significant number of mergers. 

We now turn to compare our models to the avail- 
able observational constraints. Battye & Browne (2009) and 
Browne & Battye (2010) find a higher degree of alignment 
between radio and optical emission for radio-quiet BHs 
hosted in bluer, less concentrated galaxies and a lower de- 
gree of alignment for redder, more concentrated galaxies 
and radio-loud objects. Though, they do not find any spe- 
cific dependence on stellar mass as we do (BH mass lin- 
early correlates with stellar mass), their result is in broad 
agreement with our finding on spin orientation with gas 
richness. Orban de Xivry et al. (2011) suggest that Nar- 
row Line Seyfert Is are fed through secular evolution and 
suggest how this should imply coherency in the angular 
momentum of the spin and of the disc. On the other 
hand, the degree of alignment is very low to inexistent for 
Seyfert galaxies (Kinney et al. 2000; Middelberg et al. 2004; 
Verdoes Kleijn & de Zeeuw 2005; Gallimore et al. 2006), 
which is in contradiction with our results. However, the ori- 
entation of radio jets is not necessarily aligned with BH spins 
because the propagation of jets from the accretion disc up to 
galactic scale can be strongly perturbed by the multiphase 
structure of the interstellar medium. As a consequence, gas- 
rich galaxies, which are more sensitive to gravitational in- 
stabilities with a clumpy cold gas structure, can deflect sig- 
nificantly the propagation of jets. Additionally, most well- 
studied sources with the Very Long Based Array (VLBA) 
are not disc galaxies as they are typically weak sources in 
radio, however, jet directions based on VLBA observations 
(Kellermann et al. 2004) find kinks and bends on a variety of 
scales, as well as expected systematic effects such as preces- 
sion, making it hard for available data to provide a reliable 
indication of the jet direction at the launching point, and its 
relationship with the axes of rotation of their hosts. 

Finally, we measured the distribution of BH spins as a 
function of the gas fraction of their host galaxy (Fig. 11, 
bottom panel). Spins weakly correlate with gas fractions, 
gas-rich galaxies host BHs with maximum spin and gas-poor 
galaxies host BHs with spin values around a = 0.95 slightly 
below the maximum. The reason is that gas-poor galaxies 



can be galaxies where the merger fraction in the BH mass 
build-up is any between and 1. Thus, if the galaxy has 
experienced zero mergers and the gas content is yet close 
to zero, the spin amplitude can still be close to maximum, 
as no event intervened to change it. As galaxy mergers can 
trigger strong quasar activity (Di Matteo et al. 2005), gas 
can be efficiently blown away from the central BH activity, 
and thus there is strong probability merging galaxies end up 
into gas-poor remnants. 



5 CONCLUSIONS AND DISCUSSION 

Using hydrodynamical cosmological simulations with BH 
growth and AGN feedback, we post-processed the evolution 
of BH spins with analytical prescriptions for changing spins 
through accretion of gas and BH coalescence. Our findings 
are the following: 



• Coherent gas accretion is the main driver of BH growth 
for midsize mass BHs 10 6 M < Mbh < 10 8 M Q . These BHs 
grow mostly from smoothly capturing gas from the disc of 
their host halo. Contribution from binary coalescences to 
the final BH mass becomes significant for the most massive 
BHs, A/bh > 10 8 Mq. As a consequence, midsize mass BHs 
spin to their maximum \a\ = 0.998, and the most massive 
BHs have moderate spin values \a\ ~ 0.7 — 1. 

• The ratio of accretion rate over Eddington is contin- 
uously declining over cosmic time as a results of gas rar- 
efaction in galaxies, while the contribution from mergers is 
increasing. At high redshift, even the most massive BHs have 
spin close to their maximum because the influence of merg- 
ers is transitory as gas accretion is kept at high rates and 
maintains spins aligned with their gas surroundings. As a 
consequence, the radiative efficiency of the accretion pro- 
cess onto BHs is e r = 0.32 at high redshift, and lower at low 
redshift and for the most massive BHs with e r = 0.1 for BHs 
with M BH > 10 9 M Q at z = 0. 

• AGN feedback affects spin evolution indirectly, by de- 
creasing the amount of gas available for BH accretion and its 
temperature, but not directly by homogenizing the nuclear 
gas angular momentum. 

• Spin direction does not perfectly correlate with the di- 
rection of the stellar angular momentum of its host galaxy, 
particularly for the most massive galaxies and gas-poor 
galaxies. BHs in these galaxies accrete, at very low levels, 
from the turbulent hot intra-cluster gas that bears no con- 
nection to the angular momentum of the stellar component 
of the galaxies. 

• There is a weak correlation of the BH spin amplitude 
with the gas-richness of its host galaxy. 

We have made several assumptions in our model, one of 
which is that the gas angular momentum evaluated at kpc 
scales in our simulations maintains its properties, namely its 
angular momentum direction, down to the sub-parsec scales 
close to the BH. We discussed the robustness of this assump- 
tion in Appendix A, using a 10 pc resolution high-redshift 
simulation and an isolated halo, and find that the gas an- 
gular momentum within the bulge of the galaxy shows a 
significant degree of correlation with the gas angular mo- 
mentum on galactic scales. Simulations by Hopkins et al. 
(2012) do not find such strong correlations in merging or 



BH spins in cosmological simulations 13 



bar-unstable galaxies, whereas simulations by Dotti et al. 
(2010) and Maio et al. (2013) suggest that BHs are able to 
quickly align with their gas surroundings without varying 
much in direction when embedded in dense circumnuclear 
gas discs. We discuss high-resolution test simulations in Ap- 
pendix A. 

We investigated the effect of adding some degree of in- 
coherence to the accreted gas angular momentum in Ap- 
pendix B. Even a small amount of coherence (gas is accreted 
in the semi-sphere defined by the kpc-scale gas angular mo- 
mentum, consistent with the degree of misalignment sug- 
gested by Hopkins et al. 2012) is sufficient to get spin values 
consistent with those obtained for purely coherent gas accre- 
tion. On the other hand, if the accretion process is treated 
as purely random, the lowest mass BHs have the lowest spin 
values and the highest mass BHs have the largest spin val- 
ues because in this conditions binary coalescences tend to 
spin-up BHs (Berti & Volonteri 2008). 

In our model we do not treat BH growth and AGN feed- 
back consistently with the spin evolution. The value of the 
spin in fact interferes with the Eddington limit of luminosity 
that sets a maximum accretion rate on the BH. If luminosity 
is capped at the Eddington limit, i.e., x — L the accretion 
rate on the BH scales as M — 4irG MBnm p / (e r <TTc) . In our 
simulations, we set e r = 0.1, corresponding to a BH spin 
a ~ 0.7, but for a spin a = 0.998 e r = 0.32 and the BH 
growth is slower by a factor three (King & Pringle 2006). 
In turn, this affects also the AGN energy release through 
feedback (see section 2.3) and the directions of jets. There- 
fore, there is room for potential improvement of the method 
by taking these effects into account directly into the code, 
rather than in post-processing. We note that the efficiency 
of AGN feedback is of lower importance as models of AGN 
feedback are already calibrated with lower efficiencies to ac- 
count for our inability to capture for the detailed small-scale 
physics. Jet formation at low redshift could be improved by 
using the direction of spins rather than local gas angular 
momentum. As we have seen in this paper, gas angular mo- 
mentum in the centre of galaxy clusters rapidly changes its 
direction because of the high level of turbulence in the hot 
intra-cluster gas, whereas the BH spin can keep its direction 
for several hundreds of million years. This self-consistent im- 
plementation will be the subject of forthcoming work. 
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Figure Al. Angle between the gas angular momentum of the 
bulge and the angular momentum of the disc measured at 30 pc 
and 1 kpc respectively in a zoom cosmological simulation. The 
dashed horizontal line is 90 degrees above which the gas in the 
bulge counter-rotates with the gas in the galaxy. 



Thome K. S., 1974, ApJ, 191, 507 

Tremaine S. et al., 2002, ApJ, 574, 740 

Tweed D., Devriendt J., Blaizot J., Colombi S., Slyz A., 

2009, A&A, 506, 647 
Verdoes Kleijn G. A., de Zeeuw P. T., 2005, A&A, 435, 43 
Volonteri M., Ciotti L., 2012, ArXiv e-prints, 1211.6840 
Volonteri M., Madau P., Quataert E., Rees M. J., 2005, 

ApJ, 620, 69 
Volonteri M., Sikora M., Lasota J.-P., 2007, ApJ, 667, 704 
Volonteri M., Sikora M., Lasota J.-P., Merloni A., 2012, 

ArXiv e-prints 
Wyithe J. S. B., Loeb A., 2003, ApJ, 595, 614 



APPENDIX A: COHERENCE OF ANGULAR 
MOMENTUM BETWEEN PC AND KPC 
SCALES FOR A MULTIPHASE 
INTERSTELLAR MEDIUM 

In order to test and validate the assumption that, on av- 
erage, the gas angular momentum on kpc scales is an ad- 
equate estimate of the gas angular momentum on smaller 
scales (parsec scales, which are the size of the BH scale of 
influence), we analysed the zoom-in simulation performed 
in Dubois et al. (2013) of a progenitor of a massive cluster 
at z — 6 with Ax ~ 10 pc resolution with the RAMSES code. 
The halo is simulated with the same physics than in the cur- 
rent work, i.e gas cooling, star formation, SN (type II only) 
and AGN feedback, and its mass is of 5 x 10 11 Mq at z — 6. 
More detailed information can be found in the paper stated 
above. 

In that simulation, the galaxy is composed of a bulge 
(50 pc radius) and a disc (1 kpc), which is subject to gas 
fragmentation due to Toomre instabilities, thus, clumps are 
visible most of time and are captured by the central bulge 
due to a rapid migration. Hence, it is a good case for test- 
ing the spatial coherence of gas angular momentum between 
small scales, and the kpc disc scale. The procedure is the fol- 
lowing: we compute the gas angular momentum within 30 pc 
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Figure A2. Edge-on view (top) and face-on view (bottom) of 
the gas density in the isolated halo simulation at t = 850 Myr. 



and within 1 kpc, both relative to their centre of mass, and 
we measure the angle 6 between the two vectors. Fig. Al 
shows the 6 angle as a function of redshift, and it clearly 
appears that inner gas close to the bulge has some degree 
of coherence with the gas on kpc scales. The abrupt tem- 
poral changes in the alignment of both angular momenta 
at z = 9, 8.6, 8.2, 7.9, 6.7 are due to the capture of a satel- 
lite that exerts a large gravitational torque on the central 
galaxy. The mergers at z = 9 and z = 6.7 provoke the rapid 
counter-rotation of the bulge region relative to the disc but 
last for a very short duration. 

We extended this test in convergence by setting up an 
isolated halo simulation of a M vir = 10 12 M Q with NFW 
profile that has a concentration of c = 3.5, spin of A = 0.04, 
and a gas fraction of / gas = 0.15. Gas is initially in hy- 
drostatic equilibrium in its own potential and those denned 
by the DM particles composing the dark halo. DM mass is 
resolution is MoM.res = 1.7 x 10 6 M©, and the minimum 
cell size is of Ax — 9 pc. The threshold of star formation is 
no = 250 cm -3 and a minimum temperature of 1000 K pro- 
vided by cooling through metal lines. A polytropic equation 
of state kicks in above this density threshold with polytropic 



index n — 2. Supernovae feedback and AGN feedback are 
the same as in the simulation described in the body of the 
paper. Once the simulation starts, the halo loses its initial 
pressure support through gas radiative cooling and a disc 
of gas forms due to the initial rotational support. A central 
BH is inserted after 60 Myr with an initial seed BH mass of 

2 x 10 4 M Q and reaches a final mass of 7.5 x 10 6 M Q after 1 
Gyr. 

Such a high resolution provides a multiphase description 
of the interstellar medium with cold clumps of star-forming 
gas embedded in a diffuse warm and turbulent medium. 
Fig. A2 shows the galaxy both edge-on and face-on. The 
galaxy disc shows the presence of dense gas clumps, as well 
as a dynamical thick disc. The dense clumps and the disc 
are rotationally supported, albeit with angles not precisely 
aligned with plane of the galaxy, and there is a large de- 
gree of alignment with the galaxy and, thus with the central 
BH fed through this gas. Fig. A3 highlights the BH growth 
(bottom panel) and the global coherence between angular 
momentum on small and large scales, by showing the an- 
gle between the gas angular momentum within 10 pc and 
within 1 kpc, that remains close to 10° (top panel). As a 
consequence, the angular momentum of the clumps when 
they reach the central region adds up constructively to the 
BH spin. In this simulation, we followed the spin evolution 
of the central BH, and we find that the spin is constant and 
equal to the maximum value at any time with a strong level 
of alignment with its direct gas surrounding and with gas on 
the galactic scale. The BH reaches its maximum spin value 
after 37.5 Myr after its formation time (at t = 64 Myr after 
simulation startup) , which correspond to a mass increase of 

3 times the initial BH mass during that amount of time when 
the BH grows at Eddington. This simulation is also a perfect 
test-bed for the effect of AGN feedback on the angular mo- 
mentum distribution of the galaxy gas. At time t ~ 270 Myr 
a strong outflow develops from the central AGN forming a 
large-scale galactic wind, with a bolometric luminosity from 
the central BH of L bo i = e r M B HC 2 ~ 2 x 10 44 ergs" 1 . Then, 
the accretion onto the central BH proceeds through the cap- 
ture of dense clump of gas migrating from the galactic disc 
down to the bulge (Bournaud et al. 2011; Dubois et al. 2012, 
2013) triggering powerful quasar activity with luminosity in 
the range L b oi = 10 44 - 10 45 ergs" 1 (see Fig. A3). Despite 
this strong interaction between the AGN and the galaxy gas, 
the overall coherence of the gas in the disk, between 10 pc 
and 1 kpc, is not strongly perturbed. 

As a final note, we have resolution limited to 10 pc and 
it is well possible that a more resolved turbulent spectrum 
would change this result, but complementary simulations 
performed at sub-pc resolution by Maio et al. (2013) con- 
firm that, when the accreting gas possesses a definite angu- 
lar momentum being distributed in a disc, local turbulence, 
clumpy star formation, infall from recycled stellar wind ma- 
terial and stellar feedback do not induce sufficient turbulence 
to introduce a high degree of incoherence in the gas that 
accretes on the BH, even though they do not consider feed- 
back from a central AGN source, which we instead include. 
In the case of galaxy mergers or moderately bar-unstable 
discs, instead, Hopkins et al. (2012) find that accretion onto 
the central BH is much more random though they find more 
alignment than in the pure isotropic case. In their Figure 1, 
however, most of the galaxies show a strong alignment with 
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the z-axis down to ~ 1 pc with strong departures from co- 
herency when approaching 0.1 pc. This is close to their res- 
olution limit (0.1 pc) and several smoothing lengths may in 
fact be required to sample the nuclear angular momentum 
properly. 



APPENDIX B: BH SPINS FOR DIFFERENT 
DEGREES OF GAS INCOHERENCE 

We tested the effect of adding some level of incoherence in 
the measured gas angular momentum when updating the 
BH spin through gas accretion (see Dotti et al. 2013, for 
a detailed discussion on incoherence). At each coarse time 
step, we change the measured gas angular momentum with 
a random angle uniformly distributed between and 6> max , 
where # max = 7r/2 or 6> max = n. The new direction of the gas 
angular momentum is then used to update the direction of 
the new spin. 

We see that when gas accretion is completely chaotic 
(direction of gas angular momentum is totally random), it 
leads to the formation of low-mass and midsize-mass BHs 
Mbh < 10 s M Q with low spin values of \a\ = 0.2 - 0.5, 
which is consistent with results of Fanidakis et al. (2011). 
The value of spins is increasing for more massive BHs as 
they build a large fraction of their mass through mergers. 
When the accreted gas is randomised with some level of 
coherence relative to the gas angular momentum provided 
by the simulation (pointing in the same semi-sphere as the 
gas angular momentum, consistent with the degree of mis- 
alignment suggested by Hopkins et al. 2012), values of spin 
are closer to the ones predicted in our reference model. For 
Mbh < 10 8 M Q BH spins are close to unity, \a\ = 0.9-0.95, 
and more massive BHs have lower spin values. Note that, 
in this case, values of the spin are, on average, never differ 
more than 10 per cent relative to our reference model. 

We applied an extra test, where we assume that the 
gas angular momentum accreted onto the BH is the one 
provided by the simulation as in the reference model ex- 
cept when strong merger events happen. We define a strong 
merger as the coalescence of BHs with mass ratio larger than 
1 : 10. For these pairs of BHs, gas is accreted in a chaotic 
fashion (random orientations) during 50 Myr which is a typ- 
ical duration for merger-induced starbursts, and the typical 
c-folding time of a BH accreting at the Eddington rate. Af- 
ter 50 Myr, if no more mergers happen, gas accretes with 
the angular momentum measured from the simulation. The 
result of that experiment on BH spins is similar to the ref- 
erence model: values of the spins are similar to within 1 per 
cent. We find that mergers are not responsible for the ac- 
cretion of cold gas onto BHs. The vast majority of the gas 
that contributes to BH growth is accreted while galaxies are 
not merging and are smoothly acquiring gas from their host 
halo. 
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Figure A3. Top: Inclination angle between the accreted gas an- 
gular momentum onto the BH and the gas angular momentum 
in the galaxy (black), between the accreted gas angular momen- 
tum on the BH and the angular momentum in the galaxy (red), 
and between the BH spin and the accreted gas angular momen- 
tum onto the BH (blue); middle: bolomctric luminosity; bottom: 
BH mass as a function of time for an isolated disc galaxy. The 
BH follows an initial Eddington-limited growth phase followed 
by multiple accretion of dense clump of gas that triggers strong 
quasar activity. Even though gas accretion is powered by these 
compact clumps of gas at late times, their angular momentum 
show a high degree of alignment with the BH spin and with the 
global disc angular momentum. 
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Figure Bl. Median of the spin as a function of BH mass at z = 
(black solid lines) for an angular momentum with an angle respec- 
tive to the measured gas angular momentum drawn between [0, 7r] 
(top) and [0, w/2] (bottom). Dashed lines show the 20 and 80 per- 
centiles. Red lines corresponds to no random angle (result from 
Fig 4). 



